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Introduction 



Advances have been made in the study of noncritical string theory with the discovery 
of the double-scaling limit [1] in random matrix models. The double-scaling limit correlates 
the approach of the couplings in the matrix model to certain critical values with the large- N 
limit. We give herein the solution of rectangular N x M matrix models in a 'triple-scaling' 
regime — N and M are taken to infinity independently in correlation with the approach to 
critical matrix couplings. The surface interpretation of these models in certain limits has 
been described in previous work[2,3]. The interest in these models is two-fold: (a) they 
interpolate between the rather different physics of random surfaces, and branched polymers, 
and (b) they exhibit critical behaviour involving two scaling variables. We list explicitly the 
novel results/techniques contained in the present work: 

1. We find phase transitions in correspondence to singularities of mappings from R 2 to R 2 . 

2. We use recursion relations relating polynomials for different measures on the half-line. 
The significance of this lies in the fact that the derivative terms in the differential equa- 
tions that describe critical behaviour can arise from either recursion relations relating 
different values of iV or M (or linear combinations of these two parameters). 

3. The resulting scaling behavior may be governed by: (i) two decoupled ordinary differen- 
tial equations; (ii) an ordinary differential equation and a finite difference equation; or 
(iii) two coupled partial differential equations. 

4. In case (iii) above, the equations admit weak-coupling solutions in either variable, which 
are simultaneously strong-coupling expansions in the other variable. 

Since the results and techniques given in the present paper are largely unrelated to the 
contents of Ref.'s 2,3, it can be read independently of our previous work. To put it in context, 
we mention that in work with Anderson [2, 3], we solved rectangular N x M matrix models 
in two distinct double-scaling limits, with either M (limit I)' or N — M (limit 11)^ held fixed 
as N diverges. The limit I models describe dynamically branching polymers[4,3], while the 
limit II models yield two-dimensional random surfaces [6, 2]. In Ref. 7, one of us (R.C.M.) 
solved coupled vector models in a triple-scaling limit, which has some features in common 
with the rectangular matrices solved in the present work. We wish to draw attention to 
the work of G.M. Cicuta and collaborators [8], who solved rectangular matrix models in the 
planar limit several years ago*, and to the work of J. Minahan, who has solved matrix models 
for surfaces with boundaries[9,10] — some of the equations we shall obtain in the following 
also arise in his studies (see sect. 4). 

The organization of this paper is as follows: 

t This case was also investigated in Ref.'s 4,5 with M — 1. 
J Ref.'s 6,5 analyzed this case with N — M = 0. 

* We regret that we were unaware of this work when Ref.'s 2,3 were published. 
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• Sect. 1 contains some notation, the Jacobian of interest, and a list of properties of 
the orthogonal polynomials used. We also present identities involving the recursion 
coefficients, necessary to solve the triple-scaling limit. 

• In sect. 2, we give the planar limit of the string equations. We discuss the identification 
of the critical models with singularities in two-dimensional mappings. We also give a 
brief discussion of the saddle point analysis, which provides a complementary description 
of the planar limit. 

• Triple-scaling ansdtze are introduced in sect. 3. We illustrate the construction of scaling 
solutions with an explicit example. For this example, the associated two-dimensional 
mapping can be decomposed into two separate one-dimensional singularities, and the 
resulting scaling solution involves two decoupled ordinary differential equations in two 
independent variables. 

• In sect. 4, we provide examples of new scaling equations, whose general form is a single 
ordinary differential equation which includes the second scaling variable as a param- 
eter. The scaling functions satisfy finite difference equations for shifts in this second 
variable. Two of these examples are related to the simplest intrinsically two-dimensional 
singularity, the Whitney pleat [11]. 

• Sect. 5 displays a new family of multicritical models for which i) the associated singu- 
larities are again intrinsically two-dimensional, ii) the scaling equations are two coupled 
partial differential equations in the two scaling variables, and iii) the equations admit 
two qualitatively different weak-coupling expansions. 

• Some concluding remarks and conjectures are given in sect. 6. 

1. Notation and polynomial recursion relations 



We define first matrices of the following form: 

tf=( r °t o)- « 

where T is an N x M matrix with complex entries. In the following, we will assume that 
N > M without any loss of generality. These matrices have a natural action of U(iV) xU(M), 
which can be used to bring T to the form 

T _ I X M 

\ 0(AT-M)xM 

where X M is a real diagonal MxM matrix diag(x 1; x 2 , ■ ■ ■ , x M ). Analogous statements can be 
made about matrices with real or quaternionic entries, as well as circular ensembles obtained 
by exponentiating such matrices[2]. All the steps in the following can be repeated for these 
ensembles as well, albeit with some technical complications[12,13,14]. 
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We wish to study the free energy of the partition function defined as 

Z = J dHexp(-(3trV(H 2 )), 

where dH is the Haar measure, and the matrix potential is a polynomial of the form: 
V(H 2 ) = Y.p=l%H 2 P. (Note that tr H 2k+1 = 0.) Such matrix integrals lead to theo- 
ries of random surfaces, via Feynman diagram expansions [15]. The invariance of Z under 
U(iV)xU(M) leads to the following integral: 

roo M M „ 

z= nh ex p(- 2 ^ 2 ))] riN 2(iV - M)+l n R-*?) ( 3 ) 

J °°i=l i=l l<i<j<M 

where overall factors, which do not affect the critical behavior, have been dropped. From the 
Jacobian appearing in (3), it is clear that the partition function is a function of P = iV — M, 
and M, independently. 

An obvious change of variables from x to y = x 2 simplifies the analysis of (3). It is then 
clear that we need orthogonal polynomials on a half-line, and relations between polynomials 
for measures related by the insertion of powers of y. Some of these are derived in Chihara[16]. 
Let 

d/4Z] = y Z exp(-2/3y(y)), 
where one assumes the measure has support in [0, oo).t Denote orthogonal polynomials for 
such a measure by Pn\y), such that 



"d/i[/] Pg\y) Pg\y) = 5 n , m h { £ , 



where Pn\y) = y n + lower order terms. Multiplying by y yields the usual three-term recur- 
sion relation 

V PSXv) = P { nUv) + sgipSXy) + o&P^y), (4) 

where a$ = hff/hff^. There are also simple recursion relations relating polynomials 
defined for measures with different values of / 



where 



,{l)_ K> m K '_ _ P n+lW , fi x 



t The measure will be divergent at the latter end-point for several critical models, in which 
case the following identities must be regarded as analytically continued from the convergent 
cases. 



Our scaling analysis in the following will be made in terms of the recursion coefficients, 4>n 
and 6$. The coefficients appearing in (4) are easily determined by 

One can show that 9$ and 0^ appear as recursion coefficients in the formula analogous to 
(4) for orthogonal polynomials defined for x, rather than y. 

Further useful relations between the coefficients are 

These relations are independent of the matrix potential, V(y), and their derivation is el- 
ementary, so we shall not give it here. One can think of (8,9) are discrete analogues of 
'zero-curvature' conditions (i.e., they ensure that the operations of shifting the indices / and 
n commute). Their importance in the solution of these models is quite considerable. Indeed 
in section 5, we shall find a KdV-like partial differential originating from the study of these 
relations. To this equation are associated an infinite set of triple-scaling multicritical points. 
These identities become trivial in the planar limit, with I = P, n = M, and (3 — > oo keeping 
P/P and M/P finite (see sect. 3). 

Using the properties listed above, one can derive two independent 'string equations': 

2 " ^ + M " = / Ml) V V\y) P%) P%) ( 10) 

£*23> = / MCI v\y) (p%? + ■ (ii) 

The independence of these equations is obvious since / and n are completely independent 
parameters. Using the operator techniques introduced in Ref. 17, one can reexpress these 
identities as 



r2ir j\ 

J q ^yin-id^V'^n-idx)]-^ (12) 



2n + / + 1 f 2n d\ 



i f ' 

2(3 
where 

y(n - id\) = e lX + 9(l + n- id x ) + 0(n - id x ) + e~ iX 9(n - id x )4>(n - id x ) . 

The subscript I on 1/ indicates that all of the recursion coefficients carry the same superscript 

(I) (i.e., 9(n — id x ) • 1/ = 0$). Note then that these identities only involve coefficients with 
differing subscripts. These equations are also non-trivial in the planar limit. Both of these 
features contrast them from the potential independent identities, (8) and (9). 
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2. Planar limit 



We begin by considering the planar limit, in which only the contributions of surfaces 
with spherical topology survive. Explicitly in the recursion coefficient identities above, we 
set I = P and n = M, and consider the limit f3 — > oo keeping g = N//3 and q = M//3 finite. 

We assume that the recursion coefficients have a smooth limit, 9^±f ~^ 9 an d 4>M±f ~^ 4>- 
Clearly in this case, eq.'s (8,9) vanish trivially. On the other hand, eq.'s (12,13) become 
algebraic, and can be interpreted as giving the extrema of the following functional 

W(9, 0) = U(6, <f>)-g\n9-q ln0 

where 

"2tt 



U(6, 0) = 2 / V(e lX + # + + e~ lX eci>) 

JO 27T 



Explicitly, given a potential V(y) = Yl,p=l a pU p I (^P) ■, one finds 

w*) = E? E0W- fc . d4) 

p=l p *=o v > 

Recalling the relations (7) of 6 and to the usual recursion coefficients in (4), note the 
similarity of U to the functional f2 introduced in Ref. 17. 

The critical points can be identified as follows: The planar equations take the form 

g = 9d e U{e,4>) and q = 09/7(0,0) • (15) 

In the case of hermitian matrix models[17], one finds a single equation g = g(a). In that case, 

the critical values of the coupling constants are chosen such that — for i — 0, . . . , k. 
Thus the critical points are identified as points where a — a(g) is nonanalytic, or alternatively 
as points where the map g i— > a is singular. Given this point of view, it is natural to 
identify the the critical points in the present points where the two-dimensional 

map (g, q) \— > (9, 0) becomes singular. Such points occur at the vanishing of the Jacobian 
determinant [1 1] : 

dg dg 



= \ J\ 



dq dq 



(16) 



Multicritical behavior is produced by demanding the vanishing of various higher derivatives of 
g and q. In terms of the functional W, these (multi) critical points correspond to degenerate 
extrema. In the scaling limit which will be of interest below, one does not focus on the 
degenerate point exactly. Rather the singularity is resolved by also perturbing g and q near 
their critical valuesfl]. Schematically one might find 

W « $>,(0 C - $p*0 e - 9) k > - (g c - g){9 c -9)- (q c - g)(0 c - 0) (17) 
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where pi and are integers, and Aj, arbitrary constants. The variables in (17) are denoted 
with " to indicate that they can be arbitrary linear combinations of 4>, 9, g and q. Generically 
the singularity will most naturally be parameterized in terms of various linear combinations, 
rather than 0, 9, g and q directly. This analysis of the critical points is identical to that 
for coupled vector models[7]. In that case, the functional W is replaced by an effective 
potential V, which is of direct physical relevance in those models. The coefficients of the 
polynomial terms in V are all completely independent. In contrast, (14) shows that many 
of these coefficients in W are directly related. This restricts the possible singularities that 
arise. The essential two-dimensional character of the critical points leads to a much richer 
structure than appears for ensembles of square matrices, though. We do not have a complete 
classification of all possible singularities, but in the following, we will provide various explicit 
examples which illustrate very distinct scaling behaviors. Each of these examples extend to 
an infinite family of higher multicritical models, some of which will be discussed further in 
Ref. 18. 

Before proceeding with examples of triple-scaling behavior, we would like to briefly dis- 
cuss the saddle point analysis, because it provides a more physically intuitive description of 
the planar limit, and of the critical behavior. First, we mention that using the saddle point 
approximation, the planar limit of rectangular matrix models with non-vanishing a\ and a<i 
was completely solved in Ref. [8]. The following discussion applies for general potentials. 

The partition function can be written as 

= \[{dy iy f- M eM-WViVi)]) II iVi-V,) 2 ( 18 ) 
J0 i=l l<i<j<M 

with V(y) = Ylp=i a py Jj '/ / (2p) ■ The entire integrand can be expressed as exp(— (3 2 E), where 



Zn,m 
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M M M 

( 3E({y i }) = 2j2v(y i )-^- J £ l *yi-p E HW 

i=l i=l i^j;i,j=l 



Vj\ 



The last term above arises from the Vandermonde factor in the Jacobian, and can be in- 
terpreted as a repulsive electrostatic potential between different particles at yi in a one- 
dimensional system. The second term also coming from the Jacobian, represents a further 
electrostatic repulsion away from the origin, whose strength is controlled by P = N — M. 
This repulsion combines with the external potential V to balance the Coulomb repulsion 
of the particles, and confine them to some interval [a, b] at equilibrium. We will assume 
< a < 6.T Evaluating (18) for large (3 in the saddle point approximation, requires finding 
a configuration of 'eigenvalues' satisfying 

ev {Vi) n-m i £ _i_ = . (19) 



% Wvi (3 .Jrf y% - Vj 



T The following analysis is only correct for a = if the eigenvalue density vanishes there. A 
proper treatment for a nonvanishing density at the origin is most easily done in terms of the 
original variables x rather than y[19,3]. 
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Defining a density of eigenvalues, p(y), in the usual manner [20], one finds 

v ' (y) - 9 -w =q £^-/ (z)dz - <2o) 

Solving this equation for p(y) is an example of the Hilbert problem on an arc[21]. The results 
can be expressed as constraints on the function 



- u(z) \J{z- b)(z - a^j 



2z 

where u(z) = J2k=-l ^k z ■ The constants h^, a and b are fixed by demanding that F(z) — > 
1/z as \z\ — > oo, and that F(z) remains finite as \z\ — > 0. The eigenvalue density on [a, b] is 
then given by 

p(y) = —u(y)y/(b-y)(y-a) . (21) 
7iq 

By examining the equations fixing a and b, one finds that they are identical to the planar 
string equations (15) when one makes the identification 

e = h y ^b + ^) 2 and 0=i(v / 6- v / ^) 2 • (22) 

Further using this identification, the coefficients can be written as 

L-2-k P , x 2 

h k = ~ 2 E a P+fc+2 E / ^ ( 23 ) 

which is a form similar to that for U(9, 4>) in (14). 
The aim is to evaluate the planar free energy, 

Eo = qf dyp(y) [2V(y)-(g-q)ln\y\}-q 2 [ dz dyp(z) p(y) \n\z - y\ . (24) 
J a J a 

Integrating (20), one observes that 

q J dz p(z) \n\y - z\ = V(y) - ^y^ln|y| - K, 

where one can set 

K = V(a) — - - ln|a| — q / dz p(z) \n\z — a\ . 
2 J a 



Then 



Eo = q [ dyp(y) 
J a 



+ qV(a) — - - gln|a| 



Critical behaviour corresponds to a nonanalytic dependence of the free energy Eq on the 
coupling constants in V or q. For fixed couplings, this can be expressed as nonanalyticity in 
g and q. Roughly we want 

E smg ~ (g c - gV(q c - (?) • 

Such nonanalytic behavior arises at the melting phase transition, where the eigenvalue den- 
sity (21) acquires zeroes at the boundaries, a and/or b, beyond those evident in (21). Us- 
ing eq.'s (14,15,22,23), this mechanistic description is related to (16) by observing that 
| J | = 4abu(a) u{b). 



3. Triple-scaling and Split Singularities 



In this section, we discuss triple-scaling for the case of split singularities. We use the latter 
term to describe two-dimensional singularities which can be separated into two independent 
one-dimensional singularities. In double-scaling physicsfl], many relations between quantities 
associated with critical behaviour follow from the basic ansatz that 



N 
J 



N 
J 



(25) 



where t is the double-scaling variable associated with the continuum theory. So long as 
v > 0, it is reasonable to assume that the difference equations that characterize the matrix 
model at finite N are transmuted into differential equations in the variable t as (3 — > oo. In 
triple-scaling physics, as in Ref. 7, we have two variables in the planar limit that can play 
the role of N, and thus (25) is naturally extended to 



X 

J 



X 



and 



Y 



Y 

1 



(26) 



where X and Y are linearly independent combinations of M and N determined for a given 
singularity. 

In our first example, rather than imposing simply \ J\ = as in (16), we require the 
stronger condition that the individual derivatives in the Jacobian matrix all vanish (i.e., 

W = ® = Tf$ = W = TFcj) Thus the expansion of the planar equations (15) about the 
critical point will begin at second order in both variations: A0 = <fi — C and A# = 6 — 9 C . 
One finds that when any two of the derivatives are set to vanish, the remaining two vanish 
automatically. Tuning for this critical point, therefore, imposes fewer constraints on the 
matrix potential than one might expect. The minimal potential only requires that the first 
three coefficients be nonvanishing. 

There is enough freedom in the planar equations to rescale the variables at any critical 
point to set g c = [N/f3] c = 1 and 6 C = 1. We make this choice. We also set C = y 2 since 



t Critical models with the simplest tuning to produce \J\ 
solutions — see the following sections. 



have more involved scaling 
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it is non-negative — see either (22) or (6). With these choices, the matrix potential can be 
tuned so that the desired derivatives vanish by setting: 



a\ = —3 



i + y 4 

3|/ 2 - 



«2 



3^ 



«3 



One finds that q c = [M/(3] c = y 4 (3 - y 2 )/(3y 2 - 1). C remains free — different values corre- 
spond to approaching the large- N limit with rectangular matrices of different proportions. 
g c = 1 > q c > implies that 3 > y 2 > 1. Thus the critical potential is unbounded from 
below. Expanding (15) about the critical point, to leading order 



Ag 
Aq 



32/2-1 
3 



(A0 2 + 2A(f)A9 + y 2 A6 2 ) 
(A<j) 2 + 2y 2 A(j)A6 + y 2 A6 2 ) 



(27) 



3y 2 -l 

where Ag = g — g c , Aq = q — q c . Eq. (27) can be rewritten as 



yAg ± Ag 



-{y±l){yAe±A<t>) 2 , 



(28) 



3y 2 -l' 

so one has two independent quadratic singularities. 

Eq. (28) also shows that the natural parameters in which to construct the triple scaling 
solution are yg ± q and y6 ± 0. First for (26), we define t, s by 



yN + M 
yN- M 



yN + M 

P 

yN -M 



p v ~h 



-P u ~ l s 



(29) 



The same exponents are chosen in both expressions, since the independent singularities are 
the same. With these ansatze, one can Taylor expand the differences in M and P that arise 
in eq.'s (8,9,12,13) using 



9m = ~(y + - (y - l)Sd s and d P = -y5d t - ySd s 
where 5 = f3~ v . The natural scaling ansatze for the recursion coefficients are then 



(30) 



q=2 



yo 



(p+i) 

M+n 



^m+u = vi 1 -y)- e M-\ n ] E s ) • 



(31) 



q=2 



where A^ n = 5 [n((y + l)d t + (y - l)d s ) + / y{d t + d s )\ . 

We find that v — |. The potential independent equations (8,9) yield: 

d s h2(t, s) = = d t k2(t, s) . 
10 



(32) 



Then (12,13) give 

t = 3 w — i h 2 (t) --^—dtWt). 

This result, that the scaling solution is governed by two independent Painleve I equations, 
may have been anticipated given that the critical point consists of two separate quadratic 
singularities. Such one-dimensional singularities lead to the Painleve I equation for hermitian 
matrix modelsfl]. A priori though, there is no reason that the derivative terms in the 
differential scaling equations should respect this splitting of the singularity. The crucial 
result implying this separation is (32), which is solved by k 2 = &2( s ) an d h 2 = h,2(t). 

(P) 

What is the free energy of the matrix model? Z M+ltP /Z M ^ P = (M + l)h M and (6) 
imply that 

(P) _ 2 M+hP Z M _ hP ( _ 8 2 en 

a M —k exp{-o M P ) 



J M.P 



f)( P ^ ~ Z m+i,p+i2 m ,p _ pyn / a f a nF ) 

^M,P+\^M+\,P 



(33) 



Combining these results with (30) and (31), one is able to show: First, dtd s F is a constant. 
Then up to some possible analytic terms, one finds 

k 2( s ) _j o2 r _ h 2 {t) 



d 2 F = Z^' and dfF 



Thus the nonanalytic contributions to the free energy are divided into two completely sep- 
arate pieces, both with asymptotic expansions identical to that found for ordinary two- 

5 5 

dimensional gravity[l]. The relevant expansion parameters in the present case are t~? , s~2 = 

_ 5 

02 ^ yN±M _ yN±M j 2 ^ a pp earance Q f two scaling functions should not be con- 

fused with the doubling found in hermitian matrix models[22], since in the present case they 
are functions of completely independent parameters. Furthermore, both scaling equations 
are inhomogenous in the scaling variables without the appearance of any nonuniversal coef- 
ficients^]. By rescaling h 2 , &2, t and s appropriately, the results for all different values of 
0c can be reduced to 

s = k 2 {sf - \d 2 s k 2 (s) ; d 2 F = k 2 (s) 

t = h2{tf- l -d 2 t h2{t) ; d 2 t F = h 2 {t) . 

Therefore all of these models with matrices of different proportions produce identical con- 
tinuum physics. 

The natural generalization for these split singularities would be to tune the potential 
to produce separate one-dimensional singularities of independent orders. However, such 
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singularities do not arise in the triple-scaling solution of rectangular matrix models. The 
example given here does generalize to a family of split singularities of the same order, 

yAg ± Aq oc (yA9 ± A0) n (34) 

where n is a positive integer. Implicitly we include the 'topological' model with n — 1, but 
of course this does not correspond to a singular mapping. The only assumption involved 
in producing (32) was that the scaling parameters appear with the same exponents. Since 
those results were derived from the potential independent equations (8,9), they will also 
apply for the higher singularities considered in (34). Hence one will recover k<i = &2( s ) an d 
^2 — h>2(t) for all of these critical points, and the scaling behaviour will then be governed 
by two independent ordinary differential equations of the form appearing in double-scaled 
matrix models. Note that at least up to n = 4, we found that tuning to increase the power 
of the singularity in (34) by one, only requires introducing two extra coefficients in the 
minimal potential. This fact is an example of how the interrelatedness of the coefficients in 
the functional U(6,<f>) in (14) restricts the set of singularities that can be attained in these 
matrix models. 

At the beginning, we found that g c > q c > required 3 > y 2 > 1. One may ask 
what happens at the end points of this interval. With y 2 = 3, the previous analysis and 
solution remain unchanged. It is only an exceptional point since q c = [M/f3] c = 0. This 
does not mean that the number of eigenvalues is zero or even finite. From (26), one finds 
that M = /?s(s — t)/2. Therefore the number of eigenvalues diverges, but at a much slower 
rate than either N or (3. At y 2 = 1, the previous analysis fails. One sees in (28) that one 
of yAg ± Aq vanishes to the given order in the expansion of planar equations. Hence the 
nature of the singularity changes at this point, and higher order terms must be included in 
the expansion (27) for a proper analysis. In this case, the scaling solution is similar to that 
presented in the section 5, as will be discussed in Appendix A. 

Consider the saddle point solution for the critical point solved here. Using (21) and 
(23), the eigenvalue density has two extra zeroes, one at each of the two boundaries (i.e., 
u(a) = = u(b)). The saddle point analysis leading to this result is not valid in the present 
case as is evident from (22), which requires 9 C > C . As discussed above, the scaling analysis 
is valid for <p c = y 2 > 1 = 6 C . Therefore the triple scaling solution presented here corresponds 
to a more involved analytic continuation of the matrix model than considered for the saddle 
point analysis presented in the previous section. This complication does not hold for all of 
the split singularity critical points. One finds that 9 C > C for the odd powers in (34). 

4. Parametric Solutions 

In this section we exhibit some novel critical phenomena, associated with intrinsically 
two-dimensional singularities. We begin with a general analysis of the planar limit to identify 
critical points for the simplest nontrivial potential, V(y) = a\ y/2 + a2 y 2 /4. We focus on the 
scaling solution of two higher order multicritical points. The scaling analysis for the generic 
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points will be presented in the following section. These multicritical points correspond to 
deviations from the two extreme limits of vector models and square matrix models, and only 
one of the large- N parameters, M or P, diverges while the other remains finite. The recursion 
equations, presented in section 1, lead to a coupled set of differential and finite difference 
equations. The difference equations relate scaling functions at different values of the finite 
parameter. The finite parameter also appears as an extra constant in the differential scaling 
equations. Therefore we refer to the results for these critical points as parametric solutions. 
We shall show that the new methods introduced in this paper lead to results that both 
include, and extend, known results[24,2,3,4,5]. 

For the simplest nontrivial potential, the vanishing of the Jacobian determinant (16) 
becomes 

\J\ = a\ + Aaia 2 (<p c + 9 C ) + Aa\(<p 2 c + <j> c 6 c + 6%) = . (35) 
We normalize g c = 1 = 9 C and set <p c = y 2 . In this case, (35) yields 

1 + y + y 2 1 

a i = 2 ^ — r~, — ; °2 



2y + l 2y + l 

(A second solution, which yields identical results, is produced by replacing y — > —y.) Note 
that none of the individual derivatives appearing in the Jacobian matrix vanish at this critical 
point. The numerator in a\ is a positive definite function of y, and it is straightforward to 
show a\/d2 < —3/2. This implies, in particular, that sgnai = — sgna2 for critical behaviour. 
One finds that q c = y 3 (2 + y)/(2y + 1). Requiring that g c = 1 > q c > 0, one finds two 
possible ranges for y: < y < 1 and —2 < y < — 1. Only the former range satisfies 8 C > <p c , 
while only the latter corresponds to a stable potential. We define AO, A<f>, Ag, Aq as in the 
previous section, to obtain 

(36) 

A9= %TT (A0_yAe) "%TT (A0 +2AdA( t>) ■ 

This singularity has two notable features: The linear combination q = q + yg is distinguished 
since the expansion begins at quadratic order for this parameter. Therefore one has ^| = 

= t^. Eq. (36) also singles out = — y6 as a special linear combination of the recursion 
coefficients. This is because the variations of any linear combination of g and q, except for q, 
about the critical point are linear in A0. The universal behavior is independent of the choice 
of the remaining variables, so it is best made on the basis of simplifying the calculations. At 
present, we wish to focus on exceptional points where the nature of the singularity changes. 
In the previous section, we saw that such changes may occur at the end points in the allowed 
range of y. In the present example that there are two cases of interest: i) y — and ii) 
y — — 1. This fact will be verified in the analysis of the generic points in sect. 5. 

(i) y = 0: g c = l,q c = 0, 6 C = 1, 4>c = 0, a\ = 2, a 2 = -1 
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The fact that <p c = has some important implications. First from (22), the width of the 
eigenvalue distribution in the saddle point analysis has shrunk to zero. Thus the analysis 
given there is inappropriate for this particular critical point, but it suggests that the correct 
saddle point analysis is that given in [3] for vector models. In the saddle point analysis for 
those models, M/(3 — > so that the eigenvalues accumulate at a single point since their 
electrostatic repulsion is unable to keep them apart. One can verify that the critical values 
given above correspond to one of the critical points for the vector models. This would imply 
that M should remain finite in the present case. This conclusion is further supported by 
the fact that a c = 9 c (t>c = 0. Given the usual connection in (33) of a to the free energy 
for a standard scaling solution, a c = suggests that differences in M do not give rise to 
derivatives suppressed by powers of (3 in the scaling limit. 

Returning to the expansion about the critical point, one has q = q and = 0. It is 
convenient to choose 6 and P//3 = (N — M) / (5 = g — q as the remaining variables with 
which to describe the singularity. Now (36) reduces to 

A(g -q) = -2A0 - A6 2 , Aq = -2 A9A<f> (37) 

where a term of order A0 2 has been dropped from both expressions, since in the scaling 
limit they are higher order. Rewriting (37) as Aq = A# 3 + A9A(g — q), one can recognize 
this singularity as a Whitney pleat [11]. 

Let us introduce scaling variables as in (29,30,31), 

M . , < 38 > 

and 

e£ +l) = i - exp [-i/r"dt] r^fmit) , 

= _ exp [-ip-"dt] 0-*9m(t) . 

Following the previous discussion, we have not introduced m = M + n or derivatives with 
respect to s for differences in m. From (37), we know the exponents are related as \x — 2/i, 
v — 1 — 2/i and v = 1 — 3fi, so it remains to find //. For this, we turn to (8,9), and find 

= 9m- 9m-l + d tfm, ^ 
= 9m-l(fm-l ~ fm) + dt9m-l, 

provided that fi = is, which in turn implies \x = ^, and = 0. With the vanishing of this 
last exponent, one sees that as expected M//3 does not scale, but rather (38) would give 
M = s. Thus it is inappropriate to replace discrete differences in the subscripts in (39) with 
derivatives, since the derivatives would not be suppressed by negative powers of f3, and we 
would have to keep track of arbitrarily high powers of d s . 
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Applying our ansdtze to (12,13) to complete the scaling analysis, we find 

= m + 9m(f m + f m+ i) ^ 
= t + g m + #771-1 - fm ■ 

Combining these results with (40), we derive a single equation in terms of f m only, which 
takes the surprising form 

d?f m -2fl + 2tf m = 2m-l . (42) 

This is the Painleve II equation, with a constant, and it also arises in the scaling solution 
of unitary matrix models with boundaries [9]. Its solutions have been studied in detail as 
well[10]. Note that in the present case, the constant appearing on the right hand side of (42) 
is an odd positive integer, while the constant in [9,10] is an arbitrary integer corresponding 
to the number of quark flavors coupled to the models. 

The above results are valid for any m > 1. From (5), one has 0^ = 0, so g$ = 0. Thus 
at m — 1, (40) and (41) give: = g\ + dtfi and = t + g\ — ff. Therefore f\ satisfies 

t = dth + fl (43) 

the k = 2 scaling equation in the hierarchy arising from vector models [3]. This result is 
consistent with (42), which can be written as 

{dt-2f m )(dtf m + fl-t) = 2(m-l) . 

Hence at m — 1, (42) contains the second vector equation (43). 

Since M does not scale, the relation of f m to the free energy is slightly more complicated 
than in other matrix models. One has 



eM-d P F M ) ~ %™ = n %r = n e ^ 

^M.P h^) 

m=l u m —l m=l 



which leads to 

M 
m=l 

Given the asymptotic expansion arising from (42) 

I 2m — 1 i 12m(m — 1) + 5 _5 

f m = t2 1 1 * 2-... 

4 32 

one finds 

,i M' 2 | 4M 3 + M 5 

d t F M = -Mt2 + —r l + — - — t 2 + . . . 

These results agree with those for the vector models using the analysis given in [3]. One finds, 
therefore, an a priori surprising connection between solutions of (42), and the Wronskian of 
the flows of Burger's equation[3]. 
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Therefore we have confirmed that the scaling behavior of this critical point corresponds 
to a particular scaling solution of the vector models, which describe dynamically branched 
polymers. This case generalizes to the complete hierarchy of vector critical points, which are 
characterized by <p c = = q c - Without giving the full analysis, we give the solution for the 
next critical point in this series. When g c — 1, q c — 0, 9 C — 1, 0c = 0, a\ = —a^ = 3, 03 = 1, 
we find n — v — |, (x — ^, v — 0, and 

2 2 
= m + g m (f m + fmfm+1 + f m +l - 9m-l ~ 9m ~ 9m+l) 

— t — f m + fm-l9m-l + 9mfm+l + 2fm(gm + 9m-l) ■ 
Eq. (40) remains unchanged, and combined with the above results, we find 

= t 2 f m + 6(1 - 2m) fi + {At - 5f m )fi - 4f m d t f m 

+ (At + 5fi)(d t f m ) 2 - (At - 10fi)f m d 2 f m - A(d t f m fd 2 t f m 
+ 3f m (d 2 f m ) 2 + Af m d t f m dff m - 2fldff m . 
One also finds that at m — 1 

t = f? + 3fidth+d 2 h 

which is the k — 3 vector scaling equation. There is also a topological (i.e., nonsingular) 
scaling solution for g c — 1, q c — 0, 9 C — 1, <f) c — 0, a\ = 1, with // = v — ^, /i = 1, v — 0, and 
g m = — m, / m = t. In all cases, (44) determines the free energy given the scaling solutions 
for f m . 

(ii) y = -1 : g c = l,q c =l,9 c = 1, C = 1, a x = -2, a 2 = 1 

We find very similar scaling behavior at this critical point despite the fact that it describes 
very different physics. This critical point involves square matrices since (P//3) c = g c — q c = 0, 
though they may deviate from being precisely square in the scaling limit. Since 9 C — C = 0, 
(22) indicates that the critical behavior involves the approach of the inner limit of the 
eigenvalue distribution to the origin. For P = exactly, a series of scaling solutions for such 
critical points have been studied in Ref. 25. 

Considering the expansion about the critical point, one has q = q c — gc = — P//3 and 
(f> — cf> + 0. A convenient choice of variables to complete the description of the singularity 
are q and — <p. Defining A( + ) = A(9 + 0) and A^_^ = A(9 — <fi), (36) becomes 

A(g -q) = A (+) A H , Aq = 2A (+) - ^A 2 ^ (45) 
where we have dropped terms that are higher order in the scaling limit. Rewriting (45) as 

1 Q 1 

A(g — q) — ^A^ _j + Ag, one recognizes that this critical point is again characterized 

by a Whitney pleat [11]. 

We find a scaling solution similar to that for case (i). Defining 5 = (3~s, we set 

(46) 



M 
J 



l-5 2 s 
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and 

d M+n + ^M+n = 2 " 2exp[-n*0 a ] (5 2 h«\s) + 6*h®(s) + . . .) 
^S+n - ^S+n = -2exp[-n^ s ] (S^is) + 5 2 k«\s) + . . .) . 

Hence in this case, P does not scale, so no derivatives with respect to t are introduced. Eq.'s 
(12,13) produce scaling solutions for fixed /. One finds 

4" = i( S -^ +9 4')), 

4') = -1(1 + 4*f> kf - 2d a kf ) . 

The scaling function k^ satisfies 

d 2 kf } - 2kf )3 + 2skf } = 2/ + 1, (48) 

which is the Painleve II equation with a constant, once again. The scaling parameters 
associated with M and P have traded roles in (48), as compared to (42). Since M//3 scales 

in the present case, (33) allows one to determine the free energy in terms of k^ as 

d 2 s F^ = -\{s + k[ P ^ + d s k^) . (49) 

For this critical point, the potential independent results derived from (8,9) play an auxiliary 
role, relating the scaling function at different values of I 



fr(02_o JO _ J/+l)2 , o J/+1) 



Thus the free energies at different values of P will also be related for these models. 

The above results remain unchanged at I or P = 0. This corresponds to the case of 
precisely square matrices, which were studied in Ref. 25. As noted in case (i), one can write 
(48) with / = as 

(d s -2k^)(d s k^ + k^ 2 -s) = . 

Thus the conventional solution of (48) with the asymptotic expansion kf ] = ± a V2 + . . 
leads to a trivial free energy with no higher order terms. Fortunately, the 

Painleve II equation with a nonvanishing constant admits a third solution, which for I — 
is 

1 2 8 32 

yielding 

d 2 F^ = - 1 -S+ V2 _ A s -5 + 1^-8 _ . . . . 
3 2 8 16 128 

We have checked that up to rescaling variables, the universal (i.e., nonanalytic) part of the 
free energy is identical to that found for the first critical model in Ref. 25. Note that the 
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latter analysis describes the free energy in terms of a scaling function, which satisfies Painleve 
II without a constant. 

This case generalizes to a hierarchy of critical points characterized by 9 C = <p c and g c = qc- 
At P = 0, our analysis will reproduce the results of Ref. 25. There a new series of scaling 
solutions describing random surfaces, was constructed. Their asymptotic expansions remain 
governed by the conventional differential equations arising in the analysis of Hermitian matrix 
modelsfl]. For P^O, the matrices in the present models, deviate slightly from being square. 
The scaling solutions of these matrix models are the limit II solutions discussed in [2] . The 
relation between the present analysis and that given in the previous references for both these 
random surface and the previous vector critical points is under investigation [18]. 

5. Scaling Solutions with Partial Differential Equations 



We now return, for generic values of y, to (36) which gave the expansion about the 
critical points associated with the simplest nontrivial potential. Along with = {<p — y9) and 
q = (q + yg) , we choose 9 and g to give a representative parameterization of the singularity, 
(36) becomes 

Ag = — A0 - A9 2 — A<pA9 

2y + l ( 2y+l m 

Aq = -2A^A9-^-— [ (A4 ) 2 + 3y(y+l)A9 2 ) . 

Now as long as y(y +1)^0, one can define g = g — (2y + l)q/3y(y + 1) for which the 
relevant expansion is simply 

A * = -2fTi A ^ (51) 

since the remaining quadratic terms contain Acf> and hence are higher order in the scaling 
limit. Note that (51) is the complete expansion. Cubic or higher order terms are simply not 
possible since the potential V(y) is only quadratic. Choosing 9 = 9 + (2y + l)0/3y(y + 1), 
the q expansion becomes 

A ^-^( 3 ^ + 1 > A * 2 -ft^T A ^) ■ < 52 > 

Eq.'s (51,52) characterize the generic singularities. Clearly this analysis is inappropriate for 
y = and — 1, which verifies the claim in the previous section that these are higher order 
multicritical points where the form of the singularity changes. 

The quadratic form of the singularity in (52), might lead one to conjecture that Ag 
— 5^t where 5 = [3~s. The exponents are chosen here in analogy to those which characterize 
the Painleve I equation, which arises from a quadratic singularity in the standard double- 
scaling analysis. Given this choice, (51) would lead to Ag — 5 2 s. Discrete differences in 
the large-iV parameters, (3g and (3q, would then give rise to derivatives, (3~ y dt = 5dt and 
I3~ v d s = 5 3 <9 S , respectively. We will find that these choices do indeed lead to a triple-scaling 
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solution below. An essential feature of this solution is the ratio vjv — 3. We shall show 
that with this choice, a partial differential equation very similar to the full time-dependent 
KdV equation is obtained from (8,9), without any assumptions restricting V(y). Suitable 
tunings of V{y) lead to a hierarchy of additional differential equations. Each of these new 
critical points is governed by a system of two partial differential equations, one of which is 
the KdV-like equation. 



5.1 Potential Independent Relations 

In accord with (26), we begin with a general scaling ansatz for two linearly independent 
combinations of M and P 



AM + CP 

P 

BM + DP 

P 



AM + CP' 

P . 
BM + DP 

P 



-p v ~h 



(53) 



In the scaling limit, discrete differences become derivatives as follows: 

d M = -p- y Ad t - P^Bds and d P = -fi~ v Cd t - P~°Dd s 



(54) 



Now set vjv = 3 as above, and define 5 = (3 v . We also choose 6 C = 1 and <p c = as in 
previous sections, and then introduce 



9 ( p + l ) - i _ exp 

a M+n - L eX P 



-5 {An + CI) d t - 5 3 (Bn + Dl) <9 S ] ^ S q h q (t, s), 



q=2 

-5 {An + CI) d t - 5 3 {Bn + Dl) d s j 69 s ) 

q=2 



(55) 



Substituting (55) into (8,9), we find 



C 



y=± A-C- 

One may choose the plus sign here without any loss of generality. Then 

h{t,s) =y{h 2 {t,s) + g2{s)) , 

fa{t, s) =y (hs(t, s) - ^d t h 2 {t, s) + gs{s)^j , 

and finally, one solves for k± at 0(5 5 ), but consistency between (8) and (9) requires 



d s ho = ^— ^ - df 
81 AT y 2 1 



h 2 



2 C 2 y + l 
6 y 



dfh 2 



c< 



2Ty 



2 g2d t h 2 - -d s g 2 



1 

— ( 

2 



(56) 



(57) 



(5? 



where T = AD — B C. Tis the determinant of the linear transformation relating A(M/P) 
and A{P/(3) to the scaling variables, t and s, in (53), so we require that T^O. Neglecting 
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the last two terms involving g 2 , (58) is the KdV equation with s and t playing the roles of 
the time and space coordinates, respectively. We stress that this result was obtained only 
assuming a particular ratio of the exponents in (53) and (55) (i.e., vjv = 3 and q > 2), 
without restricting the matrix potential in any way. 



5.2 Potential Dependent Results for Quadratic Potential 



We return to discussing the generic critical point associated with the quadratic potential. 



For a scaling solution of (12,13), one finds v — ^, 

C 2y + 1 



92{s) 



and 



2T y 2 

Finally h 2 must satisfy the following differential equation 



(59) 



= f- 



-92 ~ 2Cg 2 h 2 - 3C 7 



2y + V 



2y + l 



, 2 C 2 y + 1 2 
h 2 - — : — d t h 2 



6 



y 



(60) 



This equation has a form similar to the Painleve I equation, again up to terms involving g 2 , 
which introduce s dependence. 

Up until this point, we have not made use of the freedom available to rescale the variables: 
s — > zus and t — > tpt; and the functions: h 2 — > 7/12 and g 2 — > r\g% to simplify our results. 
Note that these rescalings do not provide enough freedom to completely eliminate y from the 
coefficients in (58,59,60). A suitable rescaling brings these equations to the following form: 

92 = s; 



, My + 1) n 2 , 9 

* o rS"2 + 2 



2y + l 



5 s /i2 = 



+ y + 1 

2vV + 2/ + l 



ft 



2 9 
=5-2^2 - h 2 + -d t h 2 ; 



+ -o- —92®th 2 + 

y + y + 1 



3(j/ + 1) 



2Vy 2 + y + i 



Although C and T have been eliminated from the coefficients, the last two equations still 
appear rather complicated. The true simplicity of our results is made manifest by applying 
g 2 = s, and shifting h 2 = h + (2y + l)s/ \Jy 2 + y + 1, which then yields 





= d s h 



(t + s 2 ) -h 2 + -dfh 



(61) 



2sd t h 



where dfh has also been eliminated from the KdV-like equation, by using the first equation 
above. Thus (61) shows that all values of y (which indirectly determines the proportions of 
the matrices) give the same scaling behavior. Note though that the shift between h 2 and h 
introduces a separate s dependence with a coefficient involving y. 
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Given the scaling equations (61), there are some complications in determining the free 

energy. Eq. (54) shows that both &m arid dp are both proportional to dt, to leading order. 

(P) 

From (33), the expression for a M yields 

(P) 

while that for 9 X M ' yields 

^^TT)' 12 ' (63) 

The disagreement between these expressions arises because the different approximations 
introduce different nonuniversal contributions. The divergent constant explicitly included 
in (62) is a familiar term, which is regularly dropped (as we did implicitly in section 3). 
The finite difference between (62) and (63) proportional to <?2( s ) ma Y be less familiar, but 
is simply another nonuniversal contribution. Consistent expressions for dj^dpF, dj^dpF 
and d^F can be produced by taking appropriate ratios of the recursion coefficients 9 and 
a, with the result: 

Therefore up to rescalings which can be absorbed in F, we define the universal part of the 
free energy as 

F(t,s) = j dti J 1 dt 2 h(t 2 , s) (64) 

where any possible analytic dependence on t is dropped, as well as any separate s depen- 
dence. ' 

Now consider solving the scaling equations (61). The second equation is simply solved by 
h = h{t + s 2 ). Given this form, the first equation can be solved by a perturbative expansion 
identical to that constructed for the Painleve I equation, where the expansion parameter is 



T One might think that constructing dtd s F or dgF would reveal the 'correct' s dependence. 
We were unable to find consistent definitions in general. For instance, one might expect that 

dfd s F will appear to leading order in ^^(CS^ — Adp)F, but either log (0^)^/(9^)^ 

or log (a^ ^^/(^m+i)^ appear to be equally good approximations for this quantity. Un- 
fortunately the two expressions differ at the sub leading order <5 3 , while dtd s F appears at 
order <5 4 . 
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now (t + s 2 ) 5 / 2 . Thus one finds the free energy is given by 



4 , 2^ 9, . ox 567, 2x-^ 
: i5 (t + s)2+ l6 ln(t + S '- 640 (t + s ) 2 + '" 
4s/ 5 s 2 15 s 4 



(65) 



-t 2 1 H 1 o + . . . 

15 V 2 * 8 t 2 
9, M /9 s 2 9 s 4 

+ i6 ln(() + UT-32? + 

567 _s / 5s 2 35s 4 

-»°( 2T + I? + - | + 

These results can be interpreted as follows: The free energy has the conventional topological 

_ 5 

expansion in terms of t 2, but the contribution of each topology is perturbed by terms 
(s 2 /t) m = (s/t2) 2m . In the present case, these perturbations can be resummed to give an 
overall renormalization of the topological expansion parameter of the form (t + s 2 )~2. 

It is obvious from (65)that there is a dual asymptotic solution, valid when s is large and 
t is small. This dual expansion merely has £2 replaced with s 5 , with s 2 /t — > t/s 2 . In the next 
section, we shall find that the dual expansion has a rather distinct character — neverthelsss, 
the important point is the existence of a dual expansion for these critical points. 

5.3 Higher Order Singularities 

This section discusses some results for the higher order multicritical models, for which 
the KdV-like equation governs the scaling behavior. In section 3, we described how a matrix 
potential with (2n — 1) nonvanishing coefficients could be tuned to to yield a split singularity 
of the form yAg ± Aq oc (yA9 ± A<p) n . One extra nonvanishing coefficient (for a total 
of 2n) will allow for extra tuning so that yAg + Aq can be made to vanish to order n in 
the variations, A# and Atj>3 This construction leads to the singularities where the scaling 
solution of section 5.1 is relevant. 

Explicitly for the case n — 2, the minimal potential has a\, 02, 03 and G14 nonvanishing. 
We choose g c = 1 = 6 C and C = y 2 . Then setting 

ai = -2{2y 6 + 3y 5 + 3y + 2)a 4 
a 2 = Q(y 4 + y 3 + y 2 + y + l)a 4 

a 3 = -2(2y 2 + y + 2)a 4 ( 66 ) 
1 

° 4 ~ 6y 3 + 2y 2 - 2y - 1 
yields a singularity of the form 



Ag oc A(j) 2 



Aq oc 63 A0 3 + b 2 A(f> 2 A9 + b x A0A# 2 + A^ 3 



Clearly we could have also singled out yAg — Aq, but this yields the same results. 



22 



where the notation is the same as used in (50), and b{ are constants depending on y. Note 
that in this case unlike (52), one cannot choose a linear combination 6 = + x<p to eliminate 
all of the terms mixing A0 and A6 in the singularity, although the term linear in A<ft could 
be removed. For (66), one finds that q c = y 5 (6 + 2y — 2y 2 — y 3 )/(6?/ 3 + 2y 2 — 2y — 1). 
Requiring g c = 1 > q c > 0, leaves two possible ranges for y: — 1 < y < or 1 < y < £ with 
£ = -2/3 + W 1 / 3 /3 + 100V 3 /3 « 1.59867. 

Now one may investigate a scaling solution with the ansatz given in section 5.1. The 
results presented there (56,57,58), are of course unchanged. A scaling solution of (12,13) 
requires 

1 i 

0; 



v — — 



92 



i 

S2 



9'S 







t - h 2 3 + h 2 d 2 h 2 + \{d t h 2 ) 2 - ^dfh 2 



+ ^92dfh 2 - ^g 2 h 2 2 



3 y 2 + 3y+l 2 1 y(y + 4) 3 
^2 /12 - T7T ~ , i ,o 92 



5 (y + 1) 2 * L "* 10 (2/ + I) 2 ' 

where rescalings have been made to simplify these results. A further simplification comes 

i I 
from shifting h 2 = h — ^g 2 , applying g 2 = s2, and making the following definitions: 

v — 1 re 3 3x 



y + 1 ' 



a; 3 

t S2 

20 



20 



which yields 



= T - h 3 + hdlh 



(67) 



Neglecting the last term involving a, this is the second equation in the hierarchy that arises 
for scaling solutions of hermitian matrix models[l]. With the corresponding rescalings and 
the shift of h 2) (58) becomes the KdV equation 



d G h — — dj- 



-d 2 T h 
3 T 



(68) 



Note that all of the y dependence has been absorbed in the definitions of r and a, leaving 
the scaling behavior independent of the proportions of the matrices, once again. 

One can construct a perturbative solution of these scaling equations for large r : 



r, 9 7 

F = — T3 

28 

1 



14 a Ua 6 

1 + — -o- + —^7 



7 a" 



15 r I 



81 



243 



H lnr + 

18 

1 _Z 

-T 3 



1 a' 



120 



" 1 a 

3671 + 216 
25 a 

1 + 8M 



4 

T3 



s + ■ 

243 r 2 



11 t7 4 

11664 T f 



170 a' 



505 cr 



243 r | 729 



+ . . . 



+ . . . 



+ . . . 



(69) 



These results are more complicated than the previous case, (65): a conventional topologi- 

_ 7 

cal expansion can be constructed in terms of t 3, dressed at each order by perturbations 
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(s2/t3) m . The odd powers can be absorbed in an overall renormalization of the expansion 
parameter, t — > r, as seen explicitly in (69). Perturbations with even powers remain in (69), 
and can not be incorporated into such a renormalization. In part, the latter complication 
arises because eq. (68) cannot be reduced to a linear equation as in (61). Perhaps then 
one should think of (61) as the 'topological point' for the physics described by these partial 
differential equations. 

The dual expansion at this critical point (for a 3 ^> r 2 ) is 



h =0-2 



It 3 r 2 1 r 3 105 r 4 

1 + 2j"p + r r 



9 

CT2 



128 a 6 



+ 



16 

2451 
512 ' 



51 t 393 r 2 

1 r + 

10 



5403 r 3 3489 r 4 



20 



80 



9 

(72 



+ 



6 



13 
2 



1 - 



9406 r 531501 r 2 



817 J 



+ 



6536 cr 3 1634 



16 a 
741305 r 3 

9 

(72 



+ 



229444887 r 4 
104576 ^6 



+ 



(70) 

_ 3 

Again, note that it is not possible to absorb the perturbations, now in powers of ra ? , into 

i 

a renormalization of a. There is an analogous series with h ~ — era + . . . . We have given here 
the expansion of the specific heat, not the free energy, because our methods (as mentioned 
in a preceding footnote) do not enable us to unambiguously determine the derivatives of the 
free energy with respect to s. Thus, we are not able to determine from the triple-scaling 
equations, contributions which are non-analytic in s but have a r dependence of the form 
a + br. It would be a tedious but straightforward exercise to obtain these terms using the 
techniques of [26]. Note that the dual expansions (69,70)are in terms of the renormalized 
parameters r, a. Similar expansions in terms of the bare parameters t, s can be constructed, 
as in (65). 

There are instanton solutions associated with each of these weak-coupling expansions. 
Linearizing (67,68) about (69), we find solutions of the form 



e(a, r) 



1 + a 1 



a 

2 

T3 



+ . . . exp 
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a 

ro + n—2 

T3 



+ . . . 



where r : 30 - 10r 2 + r 4 = 0, r\ = -7r (r 2 - 6)/36, and a x = (r 2 + 2)/24. At a = 0, these 
solutions reduce to those found for the k = 3 Hermitian matrix model [1]. 

Similarly, linearizing about (70), we find solutions of the form 



e(cr, r) 



a s 



1 + a x — ^ + 



3 
(72 



exp 



4 i 



T 

ra + r\ — 

fj2 



+ . . . 



where now n : 12005 - 1960r 2 + 64r 4 , r = 4^(147 - 8rf)/1029, and a x = rf/49 - 1/4. 
Note that in both instances, the prefactor and the exponential are 'dressed' by series in the 
same perturbing parameters as (69,70). We have exhibited instantons for this critical point 
to display this 'dressing'. The instantons for the other critical points found in this paper do 
not differ from the double-scaling analysis instantons. 
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Analogous behaviour is found for higher multicritical points in this series. For general 
n, the singularity may be put in the form 

Ag oc A0 n 

Aq oc A9 n+1 + b 2 ApA9 n - 1 + ... 

where '. . .' indicates terms higher order in A0. In the scaling solution, g 2 = s« while h 
satisfies a nonlinear equation similar to the k — n + 1 equation in the usual hermitian 
matrix model hierarchy. This equation will yield the conventional topological expansion in 
t but perturbed by terms of the form s» /t n + 1 J , as is evident by examining the 

planar limit. Finally, these perturbations cannot be expressed as a simple renormalization 
of the topological expansion parameter. 

6. Conclusions 



It is usually believed that the character of the large- N limit is different for theories with 
dynamical variables transforming in vector-like represenations, or in the adjoint representa- 
tion. The triple-scaling analysis given in the present work shows clearly that, at least for 
d = models, this perceived difference is a reflection of the limitations of a case-by-case anal- 
ysis. We have shown explicitly that critical phenomena ranging from polymers to random 
surfaces are described by the same equations. Random surfaces are geometrically associated 
with triangulations defined by dual graphs of Feynman diagrams [15], while polymers are 
directly associated with the Feynman graphs[3]. One intriguing result is the appearance of 
the second Painleve equation in a model of polymers. For k even, one might expect that 
the higher equations should be related to the higher mKdV flows found in unitary matrix 
models — this will be discussed elsewhere [18]. 

We have endeavoured, in the present work, to provide a detailed account of the various 
subtleties that arise in the triple-scaling analysis. The methods given in this paper seem 
inappropriate for understanding the relation, if any, of the multicritical points of triple- 
scaling models to hierarchies of integrable equations e.g., a la Drin'feld-Sokolov. Why should 
one expect such a relation? Firstly, the systems of equations obtained have a quasi-one- 
dimensional character. Secondly, a detailed examination of the examples solved in this paper 
will lead to an appreciation of the miraculous cancellations that take place in order for scaling 
behaviour in both variables to occur — thus, integrability of some sort must underlie the triple- 
scaling physics. The correct mathematical framework should, presumably, encompass the 
coupled differential equation-difference equation systems we found in sect. 4. 

We found remarkable dual weak-coupling expansions for the examples in sect. 5. The 
transition from the domain of convergence of the sub-series associated with one expansion, 
to the domain of convergence of the 'dual' expansion, is of great physical interest — to under- 
stand this, note that small s in the large t expansion, (69), corresponds to strong coupling 
in the coupling constant associated with s, and vice versa in (70). Of course, it is necessary 
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to solve the problem of the surface interpretation of the continuum theories before such a 
dual relationship can be exploited effectively. 

We now suggest an interpretation of the physics of the new 'coupling constant' associated 
with the s variable. The Liouville mode is the dynamical degree of freedom of induced 
two-dimensional gravity in background gauge. It is naturally associated with the string 
coupling constant, because of the linear dilaton-like[27] coupling of the Liouville field to the 
background curvature. In light of eq. (70), it is natural to conjecture the existence of a new 
dynamical field associated with the new expansion parameter. What might the physics of 
this new field be? Based on the surface interpretations that are suggested by finite M, N 
computations [2, 3], the new dynamics that can be investigated in the triple-scaling analysis 
would appear to involve microscopic topology fluctuations. 

How might this happen? The saddle point analysis[20] computes the physics of random 
surfaces of spherical topology. Double-scaling [1] solves the physics of random surfaces to all 
orders in the topological expansion. The finer analysis provided by triple-scaling may then 
reveal more intricate sub-structure in the physics of random surfaces. 

So much for that [28]. 

Appendix A. Another Higher Order Singularity 

At the end of section 3, we mentioned that when the split singularity degenerates at 
y 2 = 1, the scaling solution is similar to that found in section 5. We will briefly outline these 

9 1 

results. With y = 1, one has g c = \ = q c = 9 C = <p c along with a\ = —3, a 2 = 3, 03 = —5. 
Defining A( + ) = A(9 + 0) and A^_) = A(9 — <fi), the singularity may be written as 

A(» -q) = ^ - ^A 2 +) A ( _ } , A(g + q) = -3A 2 +) + |a (+) A 2 _ } . (71) 

The second term in A(g — q) is actually higher order in the scaling limit, but we explicitly 
retain it for comparison purposes below. We use the scaling ansatz presented in section 5.1, 
but the series in (55) is modified so that q > 1. We also fix A = 2, B = 0, C = 1, and D — 1 
in (53), and set vjv = 2. One finds that v — |, k\ — —h\ = — (s/2)3, and k<i = ^2 where 
/i2 satisfies 

= t - Qh 1 2 h 2 - 12h 2 2 + Adf h 2 
2d s h 2 = h\dth 2 - h\d s h\ . 

After suitable rescalings and a shift of h 2 , one finds a scaling function h = h{t + s3), which 
satisfies the Painleve I equation. This single function determines the universal part of the 

free energy as in (64). So the final result is that the free energy is identical to that found for 

4 5 

pure gravity, except with a renormalized topological expansion parameter, (t + S3) 2. 

We have been unable to generalize this scaling solution to higher order singularities. 
While the higher order split singularities all degenerate at <p c = y 2 = 1, for n > 2 the 
singularities take the form 

A(g - q) oc A^ +) A ( _ } , A(g + q) oc A n {+) + bA^ 1 A 2 _ } 
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where b is some constant. These singularities all yield parametric scaling solutions for square 
matrices, as described in section 4. Eq. (71) would have the same form with n — 2, if it were 
not for the appearance of the term in A(g — q). Thus from this point of view, (71) 

seems to be an anomalous case. 
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